Optical Emission from Aspherical Supernovae 
and the Hypernova SN 1998bw 



Keiichi Maeda^, Paolo A. Mazzali^'^'"^, Ken'ichi Nomoto^'' 



O 
O 

(N 



(N 

> 

oo 
m 



in 
o 

o 



ABSTRACT 

A fully 3D Monte Carlo scheme is applied to compute optical bolometric light curves for 
aspherical (jet-like) supernova explosion models. Density and abundance distributions are taken 
from hydrodynamic explosion models, with the energy varied as a parameter to explore the 
dependence. Our models show initially a very large degree ('^ 4 depending on model parameters) 
of boosting luminosity toward the polar (z) direction relative to the equatorial (r) plane, which 
decreases as the time of peak is approached. After the peak, the factor of the luminosity boost 
remains almost constant (~ 1.2) until the supernova enters the nebular phase. This behavior is 
due mostly to the aspherical ^^Ni distribution in the earlier phase and to the disk-like inner low- 
velocity structure in the later phase. Also the aspherical models yield an earlier peak date than 
the spherical models, especially if viewed from near the z-axis. Aspherical models with ejecta 
mass ~ IOM0 are examined, and one with the kinetic energy of the expansion 2 ± 0.5 x 10^^ 
ergs and a mass of '^^Ni ^ 0.4Mq yields a light curve in agreement with the observed light curve 
of SN 1998bw (the prototypical hyper-energetic supernova). The aspherical model is also at least 
qualitatively consistent with evolution of photospheric velocities, showing large velocities near 
the z-axis, and with a late-phase nebular spectrum. The viewing angle is close to the z-axis, 
strengthening the case for the association of SN 1998bw with the gamma ray burst GRB980425. 
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INTRODUCTION 



Lately much emphasis has been placed on the 
possible role of asymmetry in supernova (SN) ex- 
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plosions, motivating a number of multi dimen- 
sional explosion simulations both for SNe la (e.g., 
Gamezo, Khokhlov, & Gran 2005; Ropke & Hille- 
brandt 2005) and for core-collapse in massive stars 
(e.g., Proga et al. 2003; Fryer & Warren 2004; 
Janka et al. 2005; Sawai, Kotake, & Yamada 2005; 
Sekiguchi & Shibata 2005). However, there has 
been surprisingly little theoretical research on the 
emission from asymmetric supernovae (see e.g., 
Maeda et al. 2006a), despite the importance of 
this observable as a tool to test the validity of ex- 
plosion models through the comparison with ob- 
servations. In this and subsequent papers, we re- 
port our fully 3D computations of light curves and 
spectra from optical (this paper) to gamma ray 
frequencies (Maeda 2006b). 
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To date, only a scries of papers by one group 
(Hoflieh 1991, 1995; Hoflich, Wheeler, & Wang 
1999) has addressed two-dimensional computa- 
tions of optical light curves for asphcrical super- 
nova models. The interesting result emerged that 
the optical luminosity is boosted toward the polar 
(z) direction in a jet-like explosion because the 
cross section of an oblate photosphere is larger 
in this direction. They examined the early-phase 
light curve (before 100 days after the explosion). 
In the present work, we compute light curves cov- 
ering more than 1 year. 

Also, we make use of fully 3D gamma ray 
transport computations for modeling optical light 
curves. To compute an optical light curve, espe- 
cially at late phases (i.e., after 100 days), it 
is very important to follow correctly the propa- 
gation of gamma rays. These are produced by 
radioactive decays and penetrate the SN ejecta, 
depositing their energies as heat. Previous works 
did not compute optical light curves starting with 
3D gamma ray transport computations. 

For application of our light curve computations, 
we chose SN Ic 1998bw as a reference. SN Ic 
1998bw is the first supernova showing a direct ob- 
servational association with a gamma ray burst 
(GRB980425: Galama et al. 1998), opening a new 
paradigm, i.e., gamma ray burst - supernova asso- 
ciation (e.g., Hjorth ct al. 2003; Kawabata ct al. 
2003; Mathcson ct al. 2003; Mazzali, ct al. 2003; 
Stanek et al. 2003; Deng, ct al. 2005) and the 
origin of both as the end product of very massive 
stars (MacFadyen & Woosley 1999; Brown et al. 
2000; Wheeler et al. 2000). Also important is that 
SN 1998bw itself was peculiar type Ic super- 
nova. It was very bright, reaching a peak V mag- 
nitude ~ —19.1 (Galama et al. 1998, but using 
different distance modulus /i = 32.76 and extinc- 
tion Ay ~ 0.1). It showed very broad absorption 
features in the early-phase spectra around max- 
imum brightness. The light curve evolution was 
much slower than the well-studied SN Ic 19941 
(e.g., Filippenko ct al. 1995), which is thought 
to be the explosion of a low mass CO star (e.g., 
Nomoto ct al. 1994). These characteristics are in- 
terpreted in the context of spherically symmetric 
ID models as consequences of a very energetic ex- 
plosion of a massive CO star (Iwamoto et al. 1998; 
Woosley, Eastman, & Schmidt 1999). Iwamoto 
et al. (1998) derived the kinetic energy of the 



explosion i?5i = ergs = 30, the ejecta 

mass Mej = 10. 4M©, the mass of the CO star 
Mco = I3.8M0, and the main sequence mass 
Afms = 4OAf0. Later the value of the explosion 
energy was updated to be £^51 = 50 by Nakamura 
et al. (2001) to obtain a better fit to the optical 
light curve. 

However, this was not the end of the story. As 
observational data at late epochs (e.g., Patat et 
al. 2001) became available, discrepancies between 
the observations and the "spherical" hypernova 
model (e.g., Iwamoto et al. 1998) were noticed 
(e.g., McKenzie & Schaefer 1999; SoUcrman et 
al. 2000; Mazzah et al. 2001). These discrep- 
ancies could not be explained simply by modify- 
ing model parameters in the same context (e.g., 
Maeda ct al. 2003a). Something new was neces- 
sary, and one promising candidate is an aspherical 
explosion (e.g., Hoflich ct al. 1999; Maeda et al. 
2002). In particular, the late-phase spectra of SN 
1998bw are well explained by an aspherical explo- 
sion (Maeda et al. 2002, 2006a). 

The optical light curve could also be affected by 
the existence of asymmetry (Hoflich et al. 1999; 
Maeda et al. 2003a). None of the spherical hy- 
drodynamic models suggested to date can repro- 
duce consistently the entire optical light curve of 
SN 1998bw, which is covered for a period of more 
than 1 year (see e.g., Nakamura et al. 2001). Fig- 
ure 1 shows that the synthetic light curve for a 
spherical model with = 50 fits only the early 
phase, while one with E^i = 10 fits only the late 
phase (see §4 for details). Only a priori param- 
eterized spherical model fits both early and late 
phases (Chugai 2000; Maeda et al. 2003a), but no 
studies to date have investigated whether possible 
asphericity in SN 1998bw can improve the models, 
in particular with regard to the late phase optical 
light curve. 

The aim of this paper is to provide theoretical 

predictions of observable signatures of asphericity 
in supernovae/hypernovae, and to verify whether 
the peculiarities of the prototypical hypernova SN 
1998bw can be explained consistently from early 
to late phases. To add consistency to our study, 
we not only examine a light curve, but also spec- 
troscopic signatures both in early and late phases. 
The contents of the present paper are the follow- 
ing. In §2, we describe a direct method to compute 
optical and gamma ray transport in asymmetric 
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supernova ejecta. Explosion models for which we 
examine optical emission are also presented. In 
§3, radiation transport effects in aspherical su- 
pcrnovac arc presented, and the differences from 
spherical models are discussed. In §4, synthetic 
light curves and their comparison to the observed 
light curve of SN 1998bw arc presented. In §5, 
the expected photospheric velocities are compared 
with those of SN 1998bw. In §6, nebular spectra 
for the present models are presented. In §7, we 
close the paper with conclusions and discussion. 
Also in an Appendix, we describe the details of 
synthetic light curves, e.g., sensitivities on various 
assumptions. 

2. METHOD AND MODELS 

2.1. Gamma Ray and Optical Transport 

We have developed a fully 3D, energy-dependent, 
and time-dependent gamma ray transport code. 
We will describe details of the gamma ray trans- 
port scheme in a subsequent paper (Maeda 2006b). 
In what follows, only a brief summary is presented. 

The code has been developed following the indi- 
vidual packet method using a Monte Carlo scheme 
as suggested by Lucy (2005). The code follows 
gamma ray transport in SN ejecta discretised in 
3D Cartesian grids and in time steps. For the 
ejecta dynamics, we assume homologous expan- 
sion, which should be a good approximation for 
SNe la/Ib/Ic. The expansion of the ejecta is taken 
into account by dealing with the Doppler shift, 
converting at every interaction a photon packet's 
energy and direction between the comoving frame 
and the rest frame. The time delay between the 
emission of gamma rays and the final point (ei- 
ther absorption or escape out of the ejecta) is fully 
taken into account. Gamma ray lines from the 
decay chains ^'^Ni ^'^Co ^^Fe and "Ni 
^^Co — > ^^Fe are included. For the interaction of 
gamma-rays with the SN ejecta, we consider pair 
production, Compton scattering (using the Klein- 
Nishina cross section), and photoelectric absorp- 
tion (using cross sections from H to Ni). 

In view of the recent investigation by Milne et 
al. (2004) showing that not all the published ID 
gamma ray transport codes give miitually consis- 
tent results, we tested our new code by computing 
gamma ray spectra based on the (spherical) SN la 
model W7 (Nomoto, Thielemann, & Yokoi 1984), 



for which many previous studies are available for 
comparison. We compared our synthetic gamma 
ray spectra at 25 and 50 days with Figure 5 of 
Milne ct al. (2004). We found excellent agree- 
ment between our results and the spectra resulting 
from most other codes, e.g. Hungerford, Fryer, & 
Warren (2003). 

Starting with the detailed gamma ray transport 
calculations, we compute optical bolometric light 
curves. We again follow the scheme presented in 
Lucy (2005) (see also Cappellaro et al. 1997). 
This is a fully time-dependent computation, tak- 
ing into acc;ount the time delay between the cre- 
ation of a photon from the deposition of a gamma 
ray and its escape from the ejecta. The trans- 
port is solved assuming a gray atmosphere, i.e., a 
frequency-independent opacity. We test the opti- 
cal transport scheme by computing the light curve 
for the simplified SN model presented in Lucy 
(2005). For the test, we use a constant (both in 
space and time) opacity, same with the one used 
in Lucy (2005). We found good agreement be- 
tween our light curve and Lucy (2005) 's Figure 2. 
For the optical opacity for aspherical supernova 
models presented in this paper, we phenomcno- 
logically use the formula given by Chugai (2000), 
i.e., K = 0.13/[1 + (id/10)^] cm^ g-i where is 
the time in day after the explosion. The formula 
was derived by computing opacity for the 13. SM© 
CO star explosion model, including the contribu- 
tion from electron scattering and bound-free tran- 
sition, similar to the one used in Nakamura et al. 
(2001). By using the phenomenological formula, 
we miss detailed opacity distribution, and there- 
fore a detailed light curve shape. Therefore, the 
detailed fit to the early phase light curve is be- 
yond the scope of the present work (apparently, 
the late phase light curve is independent from the 
opacity prescription). This is mostly in order to 
save computational time (see also Lucy 2005). 

Each packet is therefore followed by the Monte 
Carlo method consistently from its creation as a 
gamma ray to its escape from the ejecta as either 
a gamma ray or an optical photon. In this study, 
the ejecta are mapped onto a 60'^ Cartesian grid. 
In total, 5 X 10^ photon packets (unless noted) 
with equal initial energy content are used to gen- 
erate each synthetic optical light curve. When the 
packets escape from the ejecta, they are binned 
into 135 time steps, spaced logarithmically from 2 
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to 500 days, and into 10 angular zones with equal 
solid angle from 9 = 0° to 180° (here 9 is the polar 
angle from the 2-axis). 

2.2. Models 

We use Models A, C, and F of Maeda et al. 
(2002) to compute the optical emission. These 
are the explosions of a 16Mq He star (Nomoto 
& Hashimoto 1988) with explosion energy E^i = 
Ek/W^^ ergs = 10. At the beginning of the cal- 
culations, the energy is deposited in the sphere for 
which Mr = 2AMq. The asphericity of the ex- 
plosion is obtained by depositing the energy more 
toward the polar direction (z-axis). Model F is 
a spherical model, and asphericity increases from 
Model C to Model A (see Maeda et al. 2002 for 
details) . 

To obtain (at least qualitatively) a good fit to 
SN 1998bw, we examined several models by rescal- 
ing the ejecta density and velocity self-similarly, 
so that Mej oc pv^ and Ek oc M^jV-^. According 
to previous results (Iwamoto et al. 1998; Naka- 
mura et al. 2001), we always used models with 
Mej = 1O.4M0, which is roughly expected for the 
explosion of a ~ 13. SM© CO star. 

We did not independently change the mass of 
^®Ni in the ejecta before computing light curve, 
since such a change would violate self-similarity 
in the synthetic light curve in terms of Ek and 
Mej . On the contrary, we obtained an approximate 
value of ^^Ni for each model by scaling the flux 
with the original M(^^Ni) to fit the observed lumi- 
nosity. In the future we will examine light curves 
directly using values of Mej, Ek, and M(^®Ni) ob- 
tained from hydrodynamic simulations of explo- 
sions with various parameters. Fine-tuning will 
be necessary to obtain a nice light curve, which is 
time-consuming. For example, the initial mass cut 
affects the final M(^^Ni). Accordingly, the initial 
mass cut should be fine-tuned to obtain the cor- 
rect luminosity (see e.g., Maeda & Nomoto 2003b). 
Figure 2 shows the ejecta structure of Model A. 

3. RADIATION TRANSPORT IN AS- 
PHERICAL SUPERNOVAE 

Figures 3 and 4 show the last scattering points 
of photon packets reaching an observer at any di- 
rection at a given epoch for Model A with E51 = 
20. Figure 5 shows the distribution of the total 



energy content of photon packets within position 
angle 9 ^ + d9 as& function of the position angle 
of the points of last scattering. The figures were 
obtained from computational runs with 10^ pho- 
ton packets (fewer than for the computations of 
the synthetic light curves). The fluctuation seen 
in Figure 5 is probably due to Monte Carlo noise. 

The peak date of this model is ~ 15 — 20 days, 
with some variations depending on the viewing an- 
gle (§4). Figures 3 and 4 show that at 10 days only 
photons near the surface around the 2;-axis escape 
out of the ejecta. Virtually no photons are emitted 
from the equatorial region. In other words, these 
figures consist of two sector-shaped high temper- 
ature emitting regions, not even of a spheroid. At 
20 days (near the time of peak), the last scatter- 
ing points move deeper as well as diffuse to the 
equatorial region, so that they are distributed at 
all angles. The innermost oblate spheroid (or but- 
terfly) shaped region (see also Fig. 2) is still very 
optically thick and allows no photons leakage. At 
30 days the points move even deeper. 

The peak date (^ 20 days) almost coincides 
with the epoch when the distribution of last scat- 
tering points eventually covers the entire surface 
of the ejecta. This can be understood from the 
distribution of the heating source ^^Ni. Initially, 
the ejecta density is so high that gamma ray de- 
position closely follows the distribution of ^^Ni. 
The optical photons therefore start diffusing the 
ejecta from the two poles. As time goes by, the 
diffusion time scale becomes smaller and smaller. 
It becomes almost equal to the expansion time 
scale around the time of peak. Therefore, photons 
emitted near the z-axis can diffuse to the equator 
around the peak date. This is unlike the usual 
spherically symmetric view, where the peak date 
corresponds to the date when a photon emitted 
near the center can reach the surface. This char- 
acteristic behavior, i.e., the photosphere moving 
toward the equator, implies that the light curve 
(§4) and the spectroscopic characteristics (§5) be- 
fore the peak will be differently affected by the 
viewing angle at different epochs. 

Figure 5 also shows this behavior. At epochs 
approaching the peak (i.e., between 10 and 20 
days) photons diffuse to the equator. After that 
(i.e., between 20 and 30 days), the distribution 
of the photons' last scattering points does not 
evolve significantly. Note that the angular distri- 
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bution of the photons is not spherically symmetric 
(as described by the sine curve overlapping with 
the spherically symmetric model F) even after the 
peak, but is concentrated toward the pole because 
of the high temperature resulting from a larger 
abundance of ^^Ni there. 

4. OPTICAL LIGHT CURVES 

We hereafter take the distance modulus fj, = 
32.76 (using the Hubble constant Hq = 72 km 
s"^ Mpc"i and the redshift 0.0085) and the ex- 
tinction Av = 0.1 for SN 1998bw. Figure 1 
shows optical bolometric light curves for the spher- 
ical model F. The model parameters arc as fol- 
lows: (Mej/Mo, £^51, M(56Ni)/MQ) = (10.4, 10, 
0.28), (10.4, 50, 0.40). The more energetic model 
(£■51 = 50) yields an early peak date, comparable 
to the observed one, while the less energetic model 
(£^51 = 10) peaks later than the observation. Al- 
though the rising part is not well fit even by the 
model with E51 = 50 in the present study, ex- 
tending the ^^Ni distribution out toward the sur- 
face would generate the rapid rise. (Note that in 
the present study we use the original ^^Ni distri- 
bution obtained from the hydrodynamic computa- 
tion.) On the other hand, the late time light curve 
at '-^ 50 — 500 days is more consistent with the less 
energetic model. After ~ 50 days the model with 
£51 = 50 declines more rapidly than observed. A 
model with £^51 < 10, say £^51 ~ 7 (Nakamura 
et al. 2001), would give a better fit to the late 
time light curve. This behavior of synthetic light 
curves is explained by simple scaling laws govern- 
ing the peak date tpcak oc {M^-E^'^)^/'^ (Arnett 
1982) and the late time gamma ray optical depth 

oc M^-E-^ (Clocchiatti & Wheeler 1997; Maeda 
et al. 2003a). Note that similar light curves can be 
obtained by varying the mass and energy because 
of the above scaling properties. The degeneracy 
can be in principal solved by using another infor- 
mation such as velocities (§5 and 6). 

Figure 6 shows optical bolometric light curves 
obtained for the aspherical model A. The model 
parameters are as follows: (Mej/Mg, £'51, M(^'^Ni)/M0) 
= (10.4, 10, 0.31), (10.4, 20, 0.39). The viewing 
angle is near the z-axis, within 37°. Model C 
yields a light curve almost identical to Model A, 
with a slightly delayed peak date (1 or 2 days later 
than Model A for the same parameters). 



The light curve shape of the aspherical mod- 
els is different from that of the spherical model 
F (Fig. 1). For given £^51, the aspherical model 
peaks earlier than the spherical model F. This is 
due to the combined effects of (1) low densities 
along the z-axis resulting from the large isotropic 
energy and (2) extended '"'^Ni distribution along 
the same direction. Photons therefore escape eas- 
ily through the region near the 2;-axis (see also 
Figs. 3 - 5). Furthermore, the aspherical mod- 
els yield more photons before the peak date than 
model F, while the latter shows a rather sudden 
rise around the peak with no modification of ^^Ni 
distribution. In a sense, the aspherical models nat- 
urally provide mixing of ^^Ni out to the surface 
(toward the z-axis) as a consequence of the hy- 
drodynamic process (Maeda et al. 2002; Maeda 
& Nomoto 2003b), even without additional mix- 
ing processes such as Rayleigh- Taylor instabilities. 
At late phases, on the other hand, it resembles 
the spherical model with comparable or somewhat 
smaller total kinetic energy. For example, the as- 
pherical model A with £51 = 10 yields a late time 
light curve flatter than the corresponding model F 
with £51 = 10. The decline rate is likely similar 
to the spherical model with £^51 ~ 7 (Nakamura 
et al. 2001). 

As a result, the light curve is very much im- 
proved for the aspherical models as compared with 
the spherical model F (Fig. 1). The aspherical 
models reproduce better the observations of SN 
1998bw. Although we did not try to optimize the 
fit given the uncertainty and simplifications in our 
treatment of the optical opacities (§2), an aspher- 
ical model with £51 ^10 — 20 gives a nice fit. 

Figure 7 shows how the viewing angle affects 
the light curve for Model A with £51 = 20. The 
peak is earlier and brighter for an observer at a 
smaller viewing angle 6 (angle from the z-axis). 
The ratio of the (apparent) luminosity for an ob- 
server on the z-axis to that for one on the equato- 
rial plane, is close to 4 as the SN initially emerges 
around ~ 5 days after the explosion. The ra- 
tio then rapidly decreases to ^ 1.2 at the peak 
(around ~ 20 days). After the peak and up to 
^ 40 days the ratio is almost constant and ~ 1.2, 
then it becomes unity as the whole ejecta become 
optically thin. For a less energetic model, the 
phase of constant ratio before the optically thin 
phase extends to later epochs (e.g., up to ^ 50 
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days for Model A with £^51 = 10). 

To identify what causes this behavior we exam- 
ined several models (see Appendix), and conclude 
as follows. Before the peak, the boost of the lu- 
minosity toward the ^;-direction is attributed to 
the fact that photons can only diffuse out to the 
^-direction (§3), i.e., that only the region around 
the 2;-axis is emitting photons. The cross sectional 
area of this photosphere is larger for an observer 
on the z-axis than on the r-plane, yielding a larger 
luminosity for the z-axis. As the peak date is ap- 
proached, this region continues to expand toward 
the equator, making the luminosity difference for 
different directions smaller and smaller. After the 
peak, the angular distribution of the photosphere 
does not change very much any more, and the 
luminosity difference is almost constant. In this 
period, the central disk-like region absorbs 7-rays 
preferentially along the z-axis (because of the as- 
pherical ^^Ni distribution), and emits optical pho- 
tons preferentially toward the z-direction (because 
of the disk- like structure). 

It is now interesting to compare our results to 
the previous work by Hoflich et al (1999). Al- 
though the qualitative behavior, a brighter SN for 
smaller 9, is mutually consistent, the details are 
different. First, our peak date depends on 9 since 
the extended ^^Ni distribution along the z-axis al- 
lows different diffusion time scales for different 6. 
This effect was simply not included in the pre- 
vious work. Second, before the peak the behav- 
ior is different. Our model gives initially a very 
large enhancement of the luminosity toward the 
z-axis, then the effect gradually decreases toward 
the peak date. On the other hand, Hoflich et al. 
(1999) gave an almost constant factor of the boost 
until the peak is reached, and after that the factor 
decreases by about a factor of 2 in 2 weeks (their 
Figure 6). 

The reason can be seen in the difference of 
the mechanism of boosting luminosity. Although 
in our models this is done by diffusing photons 
from z-direction to the equatorial direction, in 
Hoflich et al. (1999) boosting luminosity is always 
attributed to an oblate spheroidal photosphere. 
Note that the input models are different. Their 
parameterized models deviate greatly from spheri- 
cally symmetry, so that the axis ratio of the photo- 
sphere is always (from the innermost region to the 
outermost region) ~ 2 (see Fig. 4 of Hoflich et al. 



1999). In outer regions such as f > 10, 000 km s"""^ 
(depending on £^51) the density does not largely 
deviate from spherical symmetry in our models. 
Only below this the density contour shows a disk 
like structure (Fig. 2), but with a smaller axial ra- 
tio than theirs. Finally, in our model the difference 
in luminosities for different 9 is much smaller than 
theirs after the peak. In this phase, the boost of 
the luminosity is attributed to both the aspherical 
^^Ni distribution and the disk-like density struc- 
ture (see above), while in Hoflich et al. (1999) 
it is due only to the disk-like structure. The de- 
gree of the boost is smaller for our model than 
Hoflich et al. (1999) partly because the deviation 
of the ejecta density distribution from spherical 
symmetry is smaller in our models, although the 
direct comparison is difficult because the mecha- 
nism is different. Also, the low angular resolution 
for small 9 in the present calculations (i.e., av- 
eraging the interval 0° ~ 37°) may also smooth 
the luminosity difference to some extent. This is 
due to the current choice of an equal solid angle 
for photon binning, which is most efficient in the 
Monte-Carlo scheme. In the future, we will use 
flner bins with a larger number of photon pack- 
ets. In any case, this resolution effect is probably 
small. 

Of course, which models are more realistic is a 
different question. Although our models are based 
on hydrodynamic computations (in this sense our 
models are at least more self-consistent), differ- 
ent parameters for explosion calculations may lead 
to density structure more similar to Hoflich et al. 
(1999). Actually, most of the differences may be 
understood in terms of the different progenitors 
(the ejecta mass ~ 2 and WMq, for Hoflich et al. 
(1999) and the present study, respectively). One 
may expect larger asphericity in the inner region 
for low mass cores. 

In any case, we emphasize that it is very im- 
portant to treat the ^^Ni distribution correctly 
and the gamma ray transport in 3D, as this is 
in the present model the main cause of the char- 
acteristic behavior of the light curves. Also, we 
would suspect based on our result and analysis 
that the naive guess that larger asphericity (by a 
jet powered explosion) always yields a larger lumi- 
nosity boost toward the z-axis may not be correct. 
For example, if ^^Ni-rich blobs penetrate entirely 
through the star as an extreme case of a jet pow- 
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ered explosion, then the SN could be brighter near 
the r-plane than on the 2-axis because two blobs 
are seen from the r-plane but just one from the 
z-axis. In our model the blobs are stopped in the 
star and therefore the SN is brighter on the 2;-axis. 
In any case, light curve behavior is highly sensitive 
to density and ^^Ni distribution, and therefore one 
should always be careful when deriving aspheric- 
ity (if possible) from the light curve. One should 
always check the consistency with other informa- 
tion, e.g., spectroscopic characteristics, as we pro- 
ceed to do in subsequent sections. 

5. PHOTOSPHERIC VELOCITIES 

A model can be further constrained from its 
spectroscopic characteristics. For example, the 
evolution of the photospheric velocity provides 
important information since it scales roughly as 
Vph oc y^x/Mej, while the light curve shape de- 
pends on these parameters differently (§4). 

The concept of the photosphere needs careful 
reconsideration if the ejecta are not spherically 
symmetric. Deriving the photospheric velocity 
from observed spectra, in fact, implicitly assumes 
that the ejecta are spherically symmetric. An ab- 
sorption minimum corresponds to the line-of-sight 
velocity of a slice (perpendicular to an observer) 
where the amount of absorption takes a maximum 
value. In the spherically sjonmetric case, this ve- 
locity is exactly what is called the photospheric 
velocity, but in an asymmetric geometry this may 
not be the case. Therefore, direct comparison be- 
tween a multi dimensional model and observations 
ultimately needs multi dimensional (early phase) 
spectrum synthesis computations (e.g., Kasen et 
al. 2004), which are beyond the scope of the 
present paper. 

We give a very rough estimate of the photo- 
spheric velocity by looking at the linc-of-sight ve- 
locity of the last scattering points of optical pho- 
tons emitted toward a given observer's direction. 
Figure 8 shows the photons' energy content dis- 
tribution as functions of the line-of-sight veloc- 
ity of the last scattering points and of the epoch. 
The brightest part gives a rough estimate of the 
photospheric velocity. By definition the photo- 
spheric velocity in spherically symmetric ejecta is 
the linc-of-sight velocity of the (nearly) last scat- 
tering point moving toward the observer, while 



Figure 8 shows the line-of-sight velocity of "all" 
the last scattering points from which photons are 
emitted in the direction to the observer. There- 
fore, taking into account photons emitted aside 
the line connecting the observer and the center 
of the SN ejecta, this distribution will probably 
somewhat underestimate the photospheric veloc- 
ity. Also, the velocity is apparently sensitively af- 
fected by the optical opacity, therefore uncertainty 
in the opacity leads to uncertainty in the photo- 
spheric velocity. 

Given the above caveat, Figure 8 shows that 
the spherical model F with E^i = 50 gives a 
photospheric velocity consistent with SN 1998bw. 
Model F with the smaller energy (£^51 = 10) gives 
a velocity smaller than observed. The aspherical 
model A yields different "photospheric" velocities 
for different directions as well as different lumi- 
nosities for different directions. Especially before 
the peak there is a large velocity boost toward 
the z-axis compared to the corresponding spher- 
ical model with the same energy. Indeed, Model 
A with £^51 = 20 gives a velocity for the z-axis 
before the peak similar to that of Model F with 
E51 = 50. 

Figure 8 shows that Model A with E^i = 20 
is consistent with SN 1998bw. The model A with 
£"51 = 10 gives a bit smaller velocity before ~ 15 
days. However, it should be stated that in Model 
A with £51 = 10, the maximum velocity in the 
ejecta is ^ 20,000 km s^^. Because the decreas- 
ing density near the surface is difficult to follow 
correctly in 2D/3D hydrodynamic simulations, we 
might be underestimating the density at > 20, 000 
km s~^. After the peak the photosphere moves 
deeper and this uncertainty is not a problem. If we 
judge the model photospheric velocity at the peak, 
even Model A with £51 = 10 may marginally be 
consistent with the observation. Therefore, the as- 
pherical model A with E51 = 10 ~ 20 is consistent 
with the observed photospheric velocities, assum- 
ing the observer is placed near the z-direction. 

6. NEBULAR SPECTRA 

After the ejecta become optically thin, the SN 
evolves into the nebular phase. The emission pro- 
cess is now different from the earlier photospheric 
phase. For the nebular phase a 3D spectrum syn- 
thesis method has already been developed (Kozma 
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et al. 2005; Macda ct al. 2006a). We compute the 
nebular emission at 350 days after the explosion 
for the present models. The detail of the nebular 
emission computation in 3D and its application to 
SN 1998bw are presented in Maeda et al. (2006a). 
In Maeda et al. (2006a), the mass of some ele- 
ments are modified independently to obtain a de- 
tailed fit, while in this work we do not allow this. 
Because in this paper we rescale the hydrodynamic 
models somewhat differently, wc again compute 
nebular spectra for the present models to test con- 
sistency. 

Figure 9 shows synthetic nebular spectra for 
the present models at 350 days after the explo- 
sion as compared with the observed spectrum of 
SN 1998bw at 337 days after the peak (Patat et 
al. 2001). Although no attempt has been made 
to obtain a good fit to SN 1998bw, the aspherical 
model A is at least qualitatively consistent with 
the observation. For example, no spherical models 
give the sharply peaked [01] 6300, 6363A doublet, 
the strongest line, as is observed (e.g., Sollerman 
et al. 2000; Mazzah et al. 2001). On the other 
hand, the aspherical model A with E51 = 10 ~ 20 
yields the sharply peaked [01] line, with a view- 
ing angle close to the 2:-axis as is consistent with 
the light curve (§4), with photospheric velocity 
analysis (§5), and with the detection of a GRB. 
See Maeda et al. (2002, 2006a) for details. The 
observed nebular line broadness agrees well with 
Model A with £^51 = 10, but is a bit narrower than 
the synthetic lines of Model A with E^i = 20. This 
may possibly imply a less massive progenitor for 
SN 1998bw than examined in the present study 
(Mms = 4OM0), which is further discussed in §7. 

7. CONCLUSIONS AND DISCUSSION 

The aims of this paper are (1) to examine the 
optical properties of aspherical supernovae, in- 
cluding light curves covering more than 1 year, 
early phase photospheric velocities, and late phase 
nebular spectra, and (2) to investigate if one can 
obtain a model consistent with hypernova SN 
1998bw in all these aspects. We have presented 
observational signatures expected from aspheri- 
cal hypernova models of Maeda et al. (2002) 
and closely compare them with observations of SN 
1998bw. 

We use a Monte Carlo scheme to solve the time- 



dependent radiation transport problem in homol- 
ogously expanding SN ejecta in 3D grids. It is 
demonstrated that the scheme is very suitable to 
the problem. For example, the concept of the last 
scattering points, which can easily be traced in the 
Monte Carlo scheme, turns out to be very useful 
to understand the physical reasons of some non- 
trivial behaviors in time-dependent 3D radiation 
transport (see e.g.. Appendix). 

The light c;iirves of the aspherical models show 
characteristics different from those of spherical 
models. The aspherical model A (or C) of Maeda 
et al. (2002) shows both (1) early emergence of op- 
tical emission and (2) boosted luminosity toward 
the ^;-axis before the peak. These are mainly at- 
tributed to concentration of ^^Ni along the z-axis. 
The degree of boosting on the 2;-axis relative to the 
r-axis decreases with time, from a factor of ~ 4 to 
~ 1.2. After the peak and until the nebular phase, 
the boosting is smaller than at earlier phases. In 
this phase, the ratio of the luminosity in the ^;-axis 
to that on the r-plane remains almost constant 
at ~ 1.2. The photospheric velocity is dependent 
on the viewing angle, and is higher on the ^;-axis 
than on the r-plane. Model A with E^i = 20 gives 
photospheric velocities comparable to those of the 
spherical model F with E51 = 50. Nebular spec- 
tra are also different for different asphericity and 
different viewing angle. If the aspherical model is 
viewed on the 2;-axis, it yields sharply peaked O 
and Mg lines, and broad Fe lines (see Maeda et al. 
2006a for more detail). 

All these features are consistent with the obser- 
vations of SN 1998bw. Note that these properties 
could not be explained in spherically symmetric 
models. The highly aspherical model (A or C in 
Maeda et al. (2002)) with - lOM©, M{'^'^m) 
^ O.4Af0, and E-,i = 10 ^ 20 gives a nice fit to 
all these observations. Strictly speaking, there is 
still a small difficulty, i.e., the early phase obser- 
vations favor the larger energy E^i = 20 (§4 and 
5) although the late phase observations favor the 
smaller energy E51 = 10 (§4 and 6). However, 
given the fact that we use the "hydrodynamic" 
model, the situation is much better than in the 
spherical models (i.e., E^i — 7 vs. E^i = 50; see 
Nakamura et al. 2001). For example, if Model 
A with E51 = 10 has more high velocity mate- 
rials with some ^^Ni at V > 20,000 km s'^ than 
found in the hydrodynamic simulation (see discus- 
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sion in §5), then the problems in the early phase 
(a somewhat late peak and smaller photospheric 
velocity than observed) may be resolved. This will 
lead to an increased expansion kinetic energy, and 
£"51 ~ 20 ± 5 may be a solution. Another possi- 
bility is that the progenitor may be less massive 
than examined in the present study. A less mas- 
sive star tends to be compact in inner regions and 
to be extended in outer regions, so that it may 
yield average velocities higher in the early phases 
(i.e., outer parts) and lower in the late phases (i.e., 
inner parts) than a more massive star. However, 
it is unlikely that the progenitor's mass is very 
different from Mms = 4:0Mq since then the photo- 
spheric velocity will not follow the observations. In 
any case, it is an interesting possibility and should 
be addresses in future studies. Apparently, to do 
this it is necessary to compute hydrodynamics of 
the explosions for various progenitors rather than 
scaling the mass and energy, since the different 
density distribution will be an important factor. 

The mass of "''^Ni ^ O.4Af0 is consistent with 
the result of the aspherical explosion model of 
Maeda et al. (2002), where initial mass cut is 
set at Mr = 2AMq and all the ^^Ni is produced 
by explosive nucleosynthesis in the shock waves 
propagating the progenitor star. If the initial 
mass cut (initial remnant mass) is larger, another 
mechanism such as a massive disk wind (Mac- 
Fadyen 2003) is necessary. In either case, the nu- 
cleosynthesis products are suggested to be differ- 
ent from conventional spherically symmetric mod- 
els (e.g., Nagataki 2000; Maeda et al 2002; Maeda 
& Nomoto 2003b; Prute, Surman, & McLaugh- 
lin 2004a; Pruet, Thompson, & Hoffman 2004b). 
Some elements such as Zn, Ti are suggested to 
be enhanced. This could have a very important 
consequence on the Galactic chemical evolution 
(e.g., Kobayashi et al. 2006), especially at earliest 
phases (e.g., Iwamoto ct al. 2005). These nuclc- 
osynthetic features could be examined by observ- 
ing abundance patterns of either extremely metal 
poor halo stars (e.g., Christlieb et al. 2002: Frebel 
et al. 2005) or binary systems having experienced 
a supernova explosion of the primary star (e.g., 
Podsiadlowski et al. 2002; Gonzalez Hernandez et 
al. 2005). 

A viewing angle 9 < 30° is found to be consis- 
tent with, indeed even favorable for, the observed 
optical light curve of SN 1998bw. This is consis- 



tent with the result of Hoflich et al. (1999). The 
result seems very natural given the detection of a 
gamma ray burst in association with SN 1998bw. 
This is also consistent with a somewhat " off-axis" 
jet model {9 ~ 15°; Yamazaki, Yonetoku, & Naka- 
mura 2003) suggested for SN 1998bw. Placing 
tighter constraints on the viewing angle from mod- 
els of the optical emission is unfortunately diffi- 
cult, since the optical emission comes from non- 
relativistic ejecta which do not show relativistic 
boosting. 

Comparison of our models with other su- 
pernovae, including hypernovae SN 1997ef and 
2002ap, should be interesting. Late time light 
curves of SNe 1997ef (e.g., Mazzali, Iwamoto, & 
Nomoto 2000) and 2002ap (e.g., Mazzali et al. 
2002) are also inconsistent with spherical mod- 
els (see e.g., Maeda et al. 2003a). To solve the 
problem in these SNe will need extensive survey 
of Mej, -E51, M(^^Ni), and the ejecta geometry 
(the degree of asphericity) . Qualitatively, they 
deviate from the prediction of spherically sym- 
metric model in the same manner as SN 1998bw 
did, so that asphericity is a promising candidate 
to solve the discrepancy. If it is, then aspheric- 
ity is a general property of hypernovae. Indeed, 
the asphericity is likely a general feature of core- 
collapse supernovae. There are some supernovae, 
probably of normal energy, that possibly show a 
similar light curve behavior (e.g., Clocchiatti et 
al. 1996; Clocchiatti & Wheeler 1997). A direct 
Hubble Space Telescope image of SN 1987A clearly 
shows bipolarity (Wang et al. 2002). Also, po- 
larization and spectropolarization measurements 
suggest that core-collapse supernovae are essen- 
tially asymmetric (Wang et al. 2001; Leonard, 
Filippenko, & Ardila 2001). If available, model- 
ing not only light curves, but also nebular spectra 
(e.g., Matheson et al. 2001) is very important 
to constraint the nature of the explosion (see also 
Mazzali et al. 2005) as we see in the present paper. 

Also interesting is the high energy emission 
from aspherical supernovae/hypernovae. It should 
be noted that computations of high energy emis- 
sion based on the present model give model predic- 
tions "consistent" with optical observations of hy- 
pernova SN 1998bw. This consistent modeling be- 
tween optical and high energy emissions has only 
been done for the very nearby SN 1987A (e.g., 
McCray, Shull, & Sutherland 1987; Woosley et al. 
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1987; Shibazaki, & Ebisuzaki 1988; Kumagai ct al. 
1989). Because a hypernova is potentially a very 
interesting event not only in optical but also in X 
and 7 rays because of a large amount of radioactiv- 
ity (especially ^^Ni: see e.g., Nomoto et al. 2004) 
and large kinetic energy, model predictions consis- 
tent with the optical observations should be pro- 
vided. We will present X and 7 ray emission for 
the present models in a subsequent paper (Maeda 
2006b). 
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ciety for the Promotion of Science) Research Fel- 
lowship for Young Scientists. The authors would 
like to thank Jinsong Deng, Nozomu Tominaga, 
and Masaomi Tanaka for useful discussion, and 
Ferdinando Patat for the observational data of SN 
IQQSbw. 
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A. DETAILS OF SYNTHETIC LIGHT CURVES 



In this section, wo examine several models based on Modcil A with E^i = 20. We examine how the 
following factors affect the synthetic light curves. (1) Distribution of ^^Ni. (2) Treatment of time duration 
of a photon packet spent in the ejecta. 

Taking Model A with E^i = 20, (1) we assume that all the ''''Ni is at the center, while the density 
structure is unchanged (we call this case "central '"'^Ni"). This is an approximation to the case where 
the ^^Ni is always deeper than the photosphere. (2) We assume that the diffusion time has no angular 
dependence. We further assume Lg{t) = f{9,t)Ljnean{t), where the mean luminosity Lmean(i) is assumed 
to be that of the corresponding "spherical" model with the same Mej, M(^^Ni), and E^i, and f{0,t) is the 
angular distribution of emitted photons at fixed time t (assuming that the photon diffusion time at any angle 
is equal to that of the spherical model). These assumptions in the case (2) are used in the previous work by 
Hoflich et al. (1999) (we call this case "Approximate diffusion"). 

Following three cases are examined, (a) The central ^^Ni distribution, using the "correct" angle dependent 
diffusion time fully computed with a time-dependent 3D computation ("central ^^Ni"), (b) the original ^^Ni 
distribution, using the non-angle dependent diffusion time as described above ("Approximate diffusion"), 
and (c) the central ^^Ni and non-angle dependent diffusion, which is probably most similar to the case 
examined in Hoflich et al. (1999) ("^^Ni centered + Approximate diffusion"). Figures 10 and 11 show the 
angular luminosity distribution and light curves, respectively. Figures 12 14 show the cross sectional view 
(on the Vx — Vz plane) of the last scattering points for optical photon packets. 

How the ^^Ni distribution affects the light curve appearance is seen by comparing Figures 7a and 7b (with 
original aspherical distribution) with Figures 10a and 11a (with ^^Ni centered distribution). The effects 
of treatment of diffusion time is seen by comparing Figures 7a and 7b [10a and 11a] (directly solving 3D 
transport) with Figures 10b and lib [10c and 11c] (with the approximate diffusion). 

First, we found that the aspherical distribution of ''^'Ni, extending toward the surface along the z-axis, 
is important to make the early emergence of optical photons before the peak. All three models show the 
evolution of the light curves around the peak much slower than the original model (§4). Note that even 
Case b using the original ^^Ni distribution for computing the angle dependence is in effect similar to the 
^^Ni centered models (a and c), because of the assumption I/mean = -^spherical where the spherical model 
shows more or less centered ^^Ni distribution. Also interesting is that the models with the original ^^Ni 
distribution (original and Case b) show the boost of luminosity toward the z-axis before the peak, while the 
ones with the centered ^^Ni distribution (Cases a and c) do not. Note that while Case a apparently shows a 
luminosity boost before the peak in Figure 10a, this effect disappears well before the peak so that the effect 
is not seen in Figure 11a. A close view of the last scattering points (Figures 4 and 12 - 14) suggests that 
before the peak the boost of the luminosity toward the z-axis can mainly be attributed to the aspherical 
^^Ni distribution and energy deposition (with minor contribution due to low densities and small diffusion 
time scale along the z-axis). 

Second, we found that the treatment of the diffusion time scale is also very important. Especially we 
point out that the assumption that imean(i) is equal to the corresponding "spherical" model is not always 
correct. It ignores the fact that the ^^Ni distribution is different for the aspherical and spherical models. For 
example, comparing the original model with Case b demonstrates that the "mean" luminosity computed for 
the original model is not similar to the corresponding spherical model. 

Finally, all the three cases as well as the original model (§4) show an almost constant factor of boosting 
luminosity after the peak, and this effect is especially evident for Case a and c ('"'^Ni centered). A close view 
of the last scattering points illustrates that for Cases a and c at these epochs the photosphere is an oblate 
or disk-like spheroid with larger cross sectional area for smaller viewing angle d as suggested by Holich et 
al (1999). However, the original model and Case b with the original aspherical ^^Ni distribution show the 
distribution of the last scattering points preferentially along the z-axis. If the ^^Ni is distributed along the 
z-axis, 7-rays are absorbed preferentially along the z-axis, making temperature higher toward the z-axis. 
Because this high temperature spots are on the disk-like structure, they emit more toward the z-axis than 
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the r-axis. 

From these test calculations, we conclude that the behavior of the synthetic light curve of our aspherical 
models can be understood as follows. First, the angular dependence can be divided into two phases, i.e., 
before and after the peak. Before the peak, boosting luminosity toward the z-axis is mainly attributed 
to the aspherical ^^Ni distribution. After the peak, it is caused by combined effects of the aspherical 
^^Ni distribution and the disk-like inner structure. The fact that all the three models give light curves 
different from the original model A, especially in reproducing the rapid rise before the peak, demonstrates 
the importance of direct time dependent computations of light curves without crude approximation on the 
transport processes. 
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Fig. 1. — Light Curves for the spherical model F (Maeda et al. 2002) with the kinetic energy of the expansion 
^51 = Ek/IQ^^ = 50 (soHd) and 10 (dotted). The synthetic curves are obtained with 10^ photon packets. 
The bolometric UVOIR hght curve of SN 1998bw is taken from Patat et al. (2001), with the distance 
modulus jj, = 32.76 and the extinction Ay = 0.1. The left panel (a) shows the light curve up to 500 days, 
while the right panel (b) shows that up to only 50 days. 



15 



20000 



10000 



-10000 



-20000 



i 



20000 

-10.5 

-11 10000 



-11 -s-^ 

-12 1, 



-10000 



-20000 




0.5 



-20000 -10000 10000 20000 
Vx [cm s'^] 



-20000 -10000 10000 20000 



Vx [cm s 



Fig. 2. — Model A (Maeda et al. 2002) with E51 = 10. The left panel shows density distribution in a 
logarithmic scale (log(p [g cm"^])) at 10 days after the explosion. The ejecta is already in a homologous 
expansion phase, so that the distribution is shown in the velocity space. The right panel shows mass fractions 
of ^^Ni in a linear scale. 
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Fig. 3. — Last scattering points of optical photon packets escaping out of the ejecta in Model A with E51 = 20 

at 10 (left) and 20 (right) days after the explosion. The points are shown in the velocity space {Vx,Vy,Vz 
[cm s~^]) and the vertical axis is the Vz axis. The total number of optical photons packets is 10^. 
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Fig. 4. — Last scattering points of optical photon packets escaping out of the ejecta in Model A with £'51 = 20 
shown in the Vx — Vz plane (i.e., the slice at Vy = 0). The vertical and horizontal axes are respectively the 
Vz and the Vx axes. The last scattering points with —2 < V^/lO^cm s~^ < 2 are shown. The distribution is 
shown for 10, 20, 30, and 40 days after the explosion. 
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Fig. 5. — Angular distribution of the optical photons' total energy content. The energy content of the optical 
photons (arbitrary scaled) is binned according to the polar angle 9 of the last scattering point within the 
equal angle 4.5° each. Therefore, any spherically symmetric models should ideally give a sine curve (black). 
The distribution is shown for Model A at 10 (green), 20 (blue), and 30 days (red) after the explosion and 
for Model F (gray). The number of photons used in the computation is 10^. 
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Fig. 6. — Light Curves for the aspherical model A (Maeda et al. 2002) with the kinetic energy of the 
expansion E51 = = 20 (sohd) and 10 (dotted). The observer's direction is along the z-axis within 

9 < 37°. The synthetic curves are obtained with 5 x 10^ photon packets. The left panel (a) shows the light 
curve up to 500 days, while the right panel (b) shows that up to only 50 days. 
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Fig. 7. — Dependence of the synthetic light curve on the viewing angle 6 in Model A with = 20. The left 

panel (a) shows the synthetic light curves for different observer's directions (solid for 0° < 9 < 37°, dashed 
for 53° < 6* < 66°, and dotted for 79° < 6* < 90°). The synthetic curves are obtained with 5 x 10*" photon 
packets. The right panel (b) shows the luminosities at different normalized by the mean luminosity imean- 
Here imean is the luminosity averaged over all the solid angles (from 0° to 180°). 
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Fig. 8. — Time evolution (horizontal axes) of distribution of photon's energy contents as a function of line- 
of-sight velocities (vertical axes) at the last scattering points. The bright regions give a rough estimate of 
photospheric velocities as a function of time. See §5 for details. The distribution is shown for (a) the spherical 
model F with E51 = 50, (b) F with E51 = 10, (c) the aspherical model A with E51 = 20 and 9 = 0°, (d) A 
with £'51 = 20 and 9 = 90°, (e) A with E51 = 10 and 9 = 0°, and (c) A with E51 = 10 and 9 = 90°. Here 
9 is the polar angle (from the z-axis) of the observer's direction. Also shown is the photospheric velocity of 
SN 1998bw as a function of time (curve; Patat 2001). 
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Fig. 9. — Synthetic nebular spectra (at 350 days after the explosion) for Model F with E51 = 50 (cyan), F 
with £^51 = 10 (green), Model A with E^i — 20 (blue), and A with E^i = 10 (red). The viewing angle is 
30° from the z-axis for Model A. Also shown is the observed spectrum of SN 1998bw at 337 days after the 
maximum brightness (black). 
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Fig. 10. Synthetic light curves normahzcd by the mean himinosity Lmcan, depending on the viewing angle 
e (sohd for 0° < 6* < 37°, dashed for 53° < < 66°, and dotted for 79° < 6* < 90°) for Model A with 
i?5i = 20. Here i,„oan is the luminosity averaged over all the solid angles (from 0° to 180°). Three cases are 
shown: The panel (a) shows the model with ^^Ni distribution (artificially) concentrated at the center ("^^Ni 
centered"). The panel (b) shows the model with the assumptions imcan = -^spherical (a spherical model 
with the same energy) and no-angle dependent diffusion time ("Approximate diffusion"; See Appendix for 
details). The panel (c) shows the model with all the assumptions in (a) and (b). The models are computed 
with 10^ photon packets. 
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Fig. 11. — Synthetic light curves for the model A with £^51 = 20 with various assumptions, i.e.. (a) "•''^Ni 
centered", (b) "Approximate diffusion", and (c) "^^Ni centered and Approximate diffusion". See the caption 
of Figure 10 for details. The masses of ^^Ni in the models are O.27M0 (Case a) and O.29M0 (Cases b and 
c) (Note, however, that in Appendix we do not attempt to derive the mass of ^^Ni). 
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Fig. 12. — Last scattering points for Model A with E51 = 20, with ^^Ni distribution (artificiaUy) concentrated 
at the center ("^®Ni centered"). Last scattering points of optical photon packets are shown in the 14 — 14 
plane (i.e., the slice at = 0). The vertical and horizontal axes are respectively the 14 and the 14 axes. 
The last scattering points with —2 < V^/lO^cm < 2 are shown. The distribution is shown for 15, 20, 30, 
and 40 days after the explosion. 
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Fig. 13. — Last scattering points for Model A with E^i = 20, with the assumptions -Lmcan = isphoricai and no- 
angle dependent diffusion time ("Approximate diffusion"; See Appendix for details). Last scattering points 
of optical photon packets arc shown in the Vx — Vz plane (i.e., the slice at Vy = 0). The vertical and horizontal 
axes are respectively the Vz and the Vx axes. The last scattering points with —2 < V^/lO^cm < 2 are 
shown. The distribution is shown for 10, 20, 30, and 40 days after the explosion. 
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Fig. 14. — Last scattering points for Model A with E^i = 20 with ^^Ni distribution (artificiaUy) concentrated 
at the center, with the assumptions .Lmean = J^sphericai and no-angle dependent diffusion time ("^^Ni centered" 
and "Approximate diffusion" ; See Appendix for details). Last scattering points of optical photon packets are 
shown in the Vx — Vz plane (i.e., the slice at Vy = 0). The vertical and horizontal axes are respectively the 
Vz and the Vx axes. The last scattering points with —2 < V^/lO^cm < 2 are shown. The distribution is 
shown for 10, 20, 30, and 40 days after the explosion. 
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